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We present a lattice QCD calculation of the B — > Dlv semileptonic decay form factors f+(q 2 ) and 
fo(q 2 ) for the entire physical q 2 range. Non-relativistic QCD (NRQCD) bottom quarks and Highly 
Improved Staggered Quark (HISQ) charm and light quarks are employed together with Nf = 2 + 1 
MILC gauge configurations. A joint fit to our lattice and BABAR experimental data allows an 
extraction of the CKM matrix element V r j, \. We also determine the phenomenologically interesting 
ratio R(D) = B(B -+ Dtv t )/B(B -+ Dlv t ) (l = e,p). We find \V cb \^ t D = 0.0402(17)(13), where 
the first error consists of the lattice simulation errors and the experimental statistical error and 
the second error is the experimental systematic error. For the branching fraction ratio we find 
R(D) = 0.300(8). 

PACS numbers: 12.38.Gc, 13.20.He 


I. INTRODUCTION 

Studies of the heavy-to-heavy semileptonic decays, 
B —i Dlv and B s -+ D s lv , lead to a wealth of inter¬ 
esting and important physics. These decays can be used, 
for example, to extract the Cabibbo-Kobayashi-Maskawa 
matrix element \V c b\, providing an independent check on 
previous determinations coming from B —> D*lv decays. 
There is currently a ~ 3 a tension between the exclu¬ 
sive \V c b\ based on B -+ D*lv decays at zero recoil and 
inclusive Vj .5 determinations [T|. A recent update [5] 
by the Fermilab Lattice and MILC collaborations finds 
\V cb \^ D * = 0.03904(49) ex . pt .(53)QGD(19)QBD, whereas 

the most accurate analysis of inclusive semileptonic de¬ 
cays [3] gives \Vcb\ind . = 0.04221(78). The current un¬ 
certainty in \V c b\ leads to the dominant error in several 
important Standard Model predictions for rare decays, 
such as B s —> n + , K -+ nvv, as well as for the CP 
violation parameter ck- Reducing this uncertainty will 
have an impact on precision Flavor Physics. 

In order to get more insight into the tension between 
inclusive and exclusive |I4t>| it is crucial to determine 
\Vcb\exd. using channels other than B -+ D*lv and also 
by considering the entire physical q 2 range, rather than 
just the zero recoil point. There has been consider¬ 
able progress on this front. A very recent paper by 
the Fermilab Lattice and MILC collaborations, using 
heavy clover bottom and charm quarks, finds \ V c b\f^ D = 
(39.6 ± 1.7 QCD+exp. ± 0.2 qed) x 10~ 3 from B -+ Dlv 
lattice form factors and BABAR data [I]. And in the 
present article we give new results on B -+ Dlv form fac¬ 
tors based on the non-relativistic QCD (NRQCD) action 
for bottom and the Highly Improved Staggered Quark 
(HISQ) action for charm quarks. We also combine our 
lattice form factor results with BABAR data to extract 

\Vcb\e^d D = 0.0402(17) (13), (1) 


where the first error comes from the lattice simulation 
errors and the statistical error from experiment, and the 
second error is the systematic error from experiment. 

Interesting physics may also reside in the ratio R(D) = 
B(B -+ Dtv t )/B(B -+ Dlvi) (l = iu or e). BABAR has 
reported [5] an excess in this ratio over Standard Model 
expectations. The r lepton is considerably heavier than 
the electron or the muon, which means that the branch¬ 
ing fraction into r, v T is sensitive to both the vector form 
factor f+(q 2 ) and the scalar form factor fo(q 2 ), while the 
latter does not contribute for decays into /+ v fl or e, v e . 
This could allow scalar contributions from new physics to 
enter just in the numerator of R(D) and thereby explain 
the apparent excess. In order to confirm or reject the 
R(D) anomaly as a true new physics effect, it is impor¬ 
tant to scrutinize the current Standard Model prediction 
for R(D). Ref. |6] gave the first unquenched lattice re¬ 
sult for R(D). Using the new form factors presented in 
this article we find, 

R(D)= 0.300(8), (2) 

the most accurate Standard Model prediction to date. 

The rest of this article is organized as follows. Sec. 
II gives details of the lattice setup for this project, intro¬ 
duces the relevant bottom-charm currents and defines the 
vector and scalar form factors f+(q 2 ) and fo(q 2 )- Sec. Ill 
introduces the two- and three-point correlators that we 
simulate and describes our correlator fits and extraction 
of form factors. In Sec. IV we explain how our results for 
lattice form factors are extrapolated to the physical, i.e. 
chiral/continuum, limit. In Sec. V we discuss our form 
factor results in the physical limit and their errors com¬ 
ing from different sources. We also extract the “slope 
parameter” p 2 for f + (q 2 ). In Sec. VI we combine our 
Standard Model theory results with experimental mea¬ 
surements of the B —> Dlv branching fraction to extract 
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TABLE I. Simulation details on three “coarse” and two “fine” 
Nf = 2 + 1 MILC ensembles. 


Set 

ri/a 

mi/m s (sea) 

Nconf 

Ntsrc 

L 3 x N t 

Cl 

2.647 

0.005/0.050 

2096 

4 

24 3 x 64 

C2 

2.618 

0.010/0.050 

2256 

2 

20 3 x 64 

C3 

2.644 

0.020/0.050 

1200 

2 

20 3 x 64 

FI 

3.699 

0.0062/0.031 

1896 

4 

28 3 x 96 

F2 

3.712 

0.0124/0.031 

1200 

4 

28 3 x 96 


TABLE II. Valence quark masses aMb for NRQCD bottom 
quarks and ami and am c for HISQ light and charm quarks. 
The last column gives Z^\am c ), the tree-level wave function 
renormalization constant for massive (charm) HISQ quarks 

HU- 


Set 

aMb 

ami 

am c 

Z'f' 1 ( am c ) 

Cl 

2.650 

0.0070 

0.6207 

1.00495618 

C2 

2.688 

0.0123 

0.6300 

1.00524023 

C3 

2.650 

0.0246 

0.6235 

1.00504054 

FI 

1.832 

0.00674 

0.4130 

1.00103879 

F2 

1.826 

0.01350 

0.4120 

1.00102902 


a new value for \V c b\- Sec. VII is devoted to the ra¬ 
tio R(D). We conclude and summarize in Sec. VIII. In 
Appendix A we provide the relevant information needed 
to reconstruct our form factors, including correlations. 
Appendix B discusses further details and checks on the 
chiral/continuum/kinematic extrapolations. And in Ap¬ 
pendix C we list the priors and prior widths used in these 
extrapolations. 


II. LATTICE SETUP AND 
NRQCD/HEAVY-HISQ CURRENTS 


To study the process B -+ Dlv, one needs to eval¬ 
uate the matrix element of the bottom-charm charged 
electroweak current between the B and the D states, 
(D\(V — A) fl \B). Only the vector current V M contributes 
to the pseudoscalar-to-pseudoscalar amplitude and the 
matrix element can be written in terms of two form fac¬ 
tors f+{q 2 ) and fo{q 2 )- These depend only on the square 
of the momentum transfered between the B and the D 
mesons, q^ = p% — p f,, 


(D(p D )\V»\B(p B )} 


U(q 2 ) 


fo(q 2 ) 


Pb+Pd- 

Ml - Ml 


.M Ml-Mb . 


2- q». 


(3) 


Intermediate stages of the analysis are simplified by work¬ 
ing with the form factors /y and f±, defined by 


{D{p D )\V»\B(p B )) = v / 2Mb K/h + </J, (4) 

with 

V>1 = Mb’ P±=Pd~ ( Pd ■ v) iA (5) 

In the B rest frame (in this article we only consider B 
mesons decaying at rest) the temporal and spatial parts 
of 0 become, 

(D\V°\B) = V^/y, (6) 

(D\V k \B) = V2MBP k D f±- (7) 

Hence, one sees that one can separately determine /y or 
f± simply by looking at either the temporal or spatial 
component of V fo The conventional form factors f+(q 2 ) 
and fo(q 2 ) can then be obtained from 

u = 7m f,+ vm iMB ~ ED)h ’ (8) 


Table I lists the three coarse (a ~ 0.12fm) and two 
fine (a ss 0.09fm) MILC Nf = 2 + 1 ensembles [7] used in 
this study, together with some further simulation details. 
These MILC configurations employ the asqtad action to 
incorporate up, down and strange sea quarks. Compared 
to our recent B -+ Kl + l~ 0 0 and B s -+ Klv HD] 
projects we have increased statistics by about a factor of 
two or more. For the valence bottom quarks we use the 
NRQCD action described, for instance, in HU. The va¬ 
lence light and charm quarks are represented by the HISQ 
action |12| . In Table II we show the values for valence 
quark masses. The NRQCD bottom quark mass aMb 
was tuned in Ref. [13j to reproduce the spin averaged 
T mass, whereas the HISQ bare mass am c was tuned to 
the p c mass (suitably modified to accommodate the lack 
of annihilation and electromagnetic contributions in our 
simulations) in p3]- The valence HISQ light quark mass 
ami w &s chosen to be close to the light asqtad quark mass 
in the sea. 


fo = [{Mb ~ E ° )h + ^ ~ M ^ )M (9) 

where Eb> is the daughter D meson energy in the B rest 
frame. We generate data for four different D meson mo¬ 
menta, ff (0,0,0), |f (1,0,0), |f (1,1,0), and |f (1,1,1). 

Our goal is to evaluate the hadronic matrix elements 
(D\V°\B) and (D\V k \B) via lattice simulations. There 
are three steps in the calculation. First, one must relate 
the continuum electroweak currents, V° and V k , to lat¬ 
tice operators written in terms of the bottom and charm 
quark fields in our lattice actions. In the second step the 
matrix elements of these lattice current operators must 
be evaluated numerically and the relevant amplitudes, 
i.e. the matrix elements between the ground state B me¬ 
son and the ground state D meson with appropriate mo¬ 
menta, must be extracted. This will give us, via Eqs. 0 
and 0, the form factors /y and f± as functions of the 
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light quark mass and the D momentum. Finally, in step 
3 these numerical results must be extrapolated to the 
physical, chiral/continuum, limit. In the next two sec¬ 
tions we describe steps 2 and 3 in turn. Here we consider 
step 1 and conclude the section with a brief overview of 
the bottom-charm currents used in our simulations and of 
how these effective theory currents are matched to those 
in continuum QCD. 

Given our NRQCD action for bottom quarks and HISQ 
action for charm, we have, through next-to-leading order 
(NLO) in 1/M and lowest order in a s , the two currents 


the B meson. For instance, we create a meson at time to 
via 

= a 3 '^2'H b (x,t 0 )cl) a (x - x)j 5 %/ji(x,to)- (14) 

x' 

For the smearing functions, (f> a (x 1 — x ), we use a 5- 
function local smearing (a = 1) or Gaussian smearings 
oc e~\ x /( 2 r o) ; normalized to one (a = 2). We then 
calculate a 2 x 2 matrix of zero momentum B meson cor¬ 
relators with all combinations of source and sink smear¬ 
ings, 


4 0) =^7^, (10) 

jjl 1] = 2M b ^ >cl ‘ Jll ' V ^ b ' ( n ) 


Here ip c is the HISQ charm quark field (in its four com¬ 
ponent “naive fermion” form) and 'F b the heavy quark 
field with the upper two components given by the two- 
component NRQCD fields and the lower two compo¬ 
nents set equal to zero. We have matched these effec¬ 
tive theory currents to V/ in full QCD at one-loop or¬ 
der through 0(a s: , -%). Details of the matching 

of NRQCD/HISQ currents are given in Ref. [15]. The 
matching is similar to that employed in recent heavy-to- 
light semileptonic decays (i.e. B —► Kl + l~ and B s —► 
Klu). However, there is a difference between matching 
of NRQCD/massless-HISQ and NRQCD/massive-HISQ 
currents. Massive HISQ fermions have a nontrivial wave 


function renormalization Z^ even at tree-level. To en¬ 
sure that matching coefficients scale as {1 + 0 (a s ) +....}, 
we factor out this tree-level rescaling at the outset. This 
means the currents in (lIo| and (1II|) get multiplied by 
,- 1/2 4 

. After this rescaling, which one sees from 


(4 o) y 


Table II is a very small effect, one has 


^’“(Mo) = (15) 

x,y 

We use Gaussian widths in lattice units of size r 0 /a = 5 
on coarse ensembles and ro/a = 7 on the fine ensembles. 
For the D meson built from HISQ charm and light quarks 
we use an interpolating operator, 

& D (x,t 0 ) = a 3 'ijj c (x,t 0 )j 5 ipi(x,t 0 ), (16) 

and construct two-point correlation function with mo¬ 
mentum p, 

C D (tM\P) = ^3 KasteJ2 eiHS ~ Vl ^D(y,t)^ D ^to)). 

a?) 

The normalisation factor Ntaste = 1/16 for four- 

component naive HISQ quarks, and Ntaste = 1/4 when 
employing the one-component version of HISQ. As ex¬ 
plained in Ref. m there are no “taste” related rescaling 
factors when a nondoubled NRQCD heavy quark propa¬ 
gator is part of the loop, such as in Cg’ a (t,t 0 ) above or 
in the three-point correlator given below. 

The three-point correlator of can be written as 


(V,) QC D = (l + a sP ,)(J^)+ (12) 

with 

j(l),s ub = j (l)- as ^j(0). (13) 

Here p M and Q, are the one-loop matching coefficients 
tabulated for p = 0 and p = k in m for several aM b 
and am c values. 


III. CORRELATORS AND FITTING 
STRATEGIES 


In order to extract (D\J fi \B) (here we use “J M ” to de- 
no te e ither the full expression for the current on the rhs 
of ( 12 ) or just the lowest order term J^), we need to cal¬ 
culate the B and D meson two-point correlators and the 
three-point correlator. We use smeared heavy-light 
bilinears with Coulomb gauge fixed lattices to represent 


CJ(t,t 0 ,T-,p) = ^Y, 


Jp-(z-x) 


x (4 >d (x, to +T) Jp(z, t ) (y, to)). 


(18) 


The setup for the three-point correlator in (18) is shown 


in Fig. 1. The B meson is created at time slice A 
current insertion at timeslice t, to < t < to + T, con¬ 
verts the 6 -quark into a c-quark. The resulting D meson 
is annihilated at timeslice to + T. We have accumulated 
simulation data for four values of T: 12, 13, 14, and 15 on 
coarse and 21, 22, 23, and 24 on fine lattices. The source 
time to is picked randomly for each gauge configuration 
in order to reduce autocorrelations. Using translational 
invariance, all data are shifted to t 0 = 0 before taking av¬ 
erages and/or doing fits. The spatial sums at the source, 
Y/g, in Eqs. (fl5|), (171 and (18) are implemented using 
U( 1) random wall sources £(arjand £(x) (see Ref. [HJ for 
discussions of random wall sources in 2 - and three-point 
correlators). 
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We fit Cg’ a (t) to the form 
JVb-1 


<#“(*) = blbt*e- E ^ mt 

k=0 

N'b-1 

+ £ b*bF(-lfe 


-E- 


(19) 


k =0 


and Cn(t;p) to 


No -1 




k—0 

iV'n -1 




k—0 


( 20 ) 

The energy E E,sm differs from the full energy, E k , be¬ 
cause the NRQCD action has the 6 -quark rest mass re¬ 
moved. For the ground state, the two are related by 


E* = M B = ±(M% , --E” m )+E 0 


B,sim 


( 21 ) 


where ’ is the spin averaged T mass used to tune 
the 6 -quark mass and suitably adjusted as explained in 
Sec. II. The values of E s - m can be found in Table I of 

bb 

m. 


By comparing ( |15[ ) with (19) and (17) with (20), and 
taking the correct relativistic normalisations for the en¬ 
ergy eigenstates \E E ) and \E k ) into account, the follow¬ 
ing relations emerge, 


and 


dr — 


^EB ’ 


( 22 ) 


(23) 


For the three-point correlator Cj(t,T;p) we use the 
following fit ansatz 


Nb — 1 No — 1 


c a j{t,T-p)= e E A 


jk 


e j e k v 1 


j=0 k—0 


Nd — 1 N' b — 1 

E E 

j—0 k—0 

N d — 1 N b — 1 


a . e - E ?- t e- E '*' Sim < T - t ')<-^ T - t '> 


jk 


( — 1 ) 


E E c 




(- 1 )‘ 


j—0 k—0 

N’ b - 1N’ d -1 

E E (24) 

j—0 k—0 



FIG. 1. Setup for three-point correlators 
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FIG. 2. aEo versus N exp for several momenta and for en¬ 
semble FI. 


The amplitudes A r - k etc. depend on the current J M and 
on the D meson momentum p. Again by comparing (18) 
and (24) and using (22) and (23), one finds 


A« k = dj 


(E?\J,\E e ) 

y/2a3EPyj2a3ES 


bT- 


(25) 


For j = k = 0, Aq 0 in (25) gives us the sought after 
hadronic matrix elements (D\J k \B), 

A’ 


W„\B) = EE j2a?E? y/2a?M B . (26) 
"0 0 0 v 


do b\ 


Our fitting strategies based on Bayesian methods have 
been developed and refined in a number of calcula¬ 
tions PH EH- We follow closely the approach used in 
our recent B s —>• I\lu studies m- Figs. 2 and 3 show fit 
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FIG. 3. aEg’ alm versus N exp for ensemble FI. 


FIG. 5. A 0 o for J 0 (p = (0,0,0)) versus N D ^ B and N' DB 
for FI. The two plots on the right are at fixed Nb,b = 4 or 
N'd,b = 4. 



ap 


FIG. 4. Dispersion relations on ensembles Cl and FI. The 
shaded region is bounded by 1 ± ca s {ap) 2 with c = 0.1. 


results for the groundstate D meson energies aE B and 
for aE^’ slm , respectively, versus the number of exponen¬ 
tials in the fit N exp (we set N exp = N D ^ B = N' D B ). Fits 
stabilize after N exp = 4. In Fig. 4 we check the ratio 
{M 2 +p 2 )/E 2 for the D meson on one coarse (Cl) and 
one fine (FI) ensemble. The shaded area is bounded by 
1 ± c a s (ap ) 2 , where the free parameter c has been set 
to 0.1. One sees that the relativistic dispersion relation 
holds within errors to better than 0.5%. 

For fixed D momentum, the combination on the rhs of 



CN (N (N <N (N 
<N 


FIG. 6. Aoo for Jo{p = (0,0,0)) versus different T combina¬ 
tions for FI. We take the result with T = (21, 22, 23, 24). 

Eq. (24) is obtained from a simultaneous fit to a 2 x 2 
matrix of B-correlators, one D correlator and numerous 
three-point correlators. The number of three-point cor¬ 
relators in the fits varies from six to sixteen as we use 
either three or four values of T, two smearings a , and ei¬ 
ther one current at zero momentum (Jo) or two currents 
at non-zero momenta (Jo and Ji). We call this type of 
fit an “individual fit”. Fig. 5 shows results for for 
J 0 {p = (0, 0,0)) versus the number of exponentials. And 
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FIG. 7. Comparison between form factor fit results when 
current matching corrections are undertaken before or after 
the fits (for ensemble FI). 


Fig. 6 shows how results depend on choices for different T 
combinations. Individual fits give stable and consistent 
results under such variations. 

We fit data after the current matching described in 
Sec. II. Specifically, we obtain simulation data for J^ 
and of Eqs. (8) and (9), reconstruct the full expres¬ 
sion on the RHS of Eq (10), and fit the resulting data. 
Alternatively, we can perform separate fits to J^ and 
J^ and then combine the results according to Eq. (10). 
We have compared the two approaches and find good 
agreement, as shown in Fig. [7] 

In order to get correlations of form factors at differ¬ 
ent q 2 ' s, one needs to do a simultaneous fit with all dif¬ 
ferent D momenta. Each individual fit alone involves 
10 to 20 correlators (depending on the combination of 
T values) and a simultaneous fit requires a four times 
larger set of correlators. We find that the simultaneous 
fits lead to unreliable results, indicating that they are 
too complicated given the accuracy of our data. Our 
fitting routines, however, allow for an alternate way to 
keep track of correlations between form factors at differ¬ 
ent q 2, s. One can do a sequence of individual fits, one 
after the other, all within a single script and always em¬ 
ploying the full covariance matrix for all the data (all D 
meson momenta). We call such fits “master fits”. These 
master fits are easier than straight simultaneous fits, but 
still highly non-trivial and time consuming. It was pos¬ 
sible to get good master fits consistent with individual 
fit results; however, these were less stable with respect 
to changes in N exp and T combinations. Hence for our 


FIG. 8. Correlations between different momenta from the 
master fit for ensemble FI. 


TABLE III. Fit results for fo(p) and /+(p). 


Set 

fo (0,0,0) 

/o(l, 0,0) 

/o(l, 1,0) 

/o(l, 1 , 1) 

Cl 

C2 

C3 

0.8810(56) 

0.8809(31) 

0.8872(23) 

0.8743(43) 

0.8716(54) 

0.8685(32) 

0.8608(38) 

0.8617(44) 

0.8592(29) 

0.8534(42) 

0.8503(50) 

0.8473(38) 

FI 

F2 

0.9034(31) 

0.9051(23) 

0.8771(42) 

0.8895(36) 

0.8643(41) 

0.8702(29) 

0.8479(56) 

0.8504(34) 


Set 

/+(1,0,0) 

/+(i,i,o) 

/+(1,1,1) 

Cl 

1.135(12) 

1.1125(57) 

1.0837(61) 

C2 

1.110(12) 

1.0809(70) 

1.0479(64) 

C3 

1.1282(71) 

1.0937(40) 

1.0569(50) 

FI 

1.1344(91) 

1.0931(59) 

1.0480(74) 

F2 

1.1461(72) 

1.0963(39) 

1.0577(45) 


final fit results we use central values and errors from in¬ 
dividual fits, and use the good master fits just to extract 
the necessary correlations. Fig. [8] shows an example of 
correlations obtained from a master fit to all the data 
from ensemble FI. Table III summarizes form factor re¬ 
sults for each ensemble and D momentum, and Table 
IV shows the fit results of B and D meson ground state 
energies. 

IV. CHIRAL, CONTINUUM AND KINEMATIC 
EXTRAPOLATION 

In this section we describe how we extrapolate the form 
factors of Table III to the continuum and chiral limits, 
and how one can get information on the form factors 
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TABLE IV. Fit results for aEo with each momentum and 

rr\B ,sim 

aE 0 


Set 

aEn( 0, 0, 0) 

aEo(l, 0,0) 

aEo(l, 1,0) 

aE n ( 1,1,1) 

Cl 

1.1388(15) 

1.1681(19) 

1.1979(17) 

1.2267(16) 

C2 

1.1577(14) 

1.1959(29) 

1.2365(28) 

1.2758(34) 

C3 

1.16355(69) 

1.2046(11) 

1.2442(12) 

1.2824(18) 

FI 

0.81409(42) 

0.84349(77) 

0.87262(83) 

0.9003(13) 

F2 

0.81999(37) 

0.85051(69) 

0.87905(66) 

0.90682(86) 

: : 

Set 




tt\B ,sim 

aE 0 

Cl 




0.4964(13) 

C2 




0.5089(14) 

C3 




0.51376(95) 

FI 




0.38190(89) 

F2 




0.38726(73) 


for the entire physical kinematic range. In the contin¬ 
uum physical theory, form factors are functions of a sin¬ 
gle kinematic variable which can be taken to be q 2 , Ed, 
(w — 1) = (g^ OI - q 2 )/[2 MbMd) or the 2 - variable defined 
in terms of q 2 as, 


/ 2\ \E+ d 2 \/^+ ^0 

z(q ) = . - & 

V f + - Q 2 + Vt+ ~to 


(27) 


Here t + = (Mg + Md) 2 and to is a free parameter which 
we set to t 0 = q^ ax = ( Mb — Md) 2 ~ 11.66 GeV 2 . A 
popular expansion in terms of 2 is the Bourrely-Caprini- 
Lellouch (BCL) parametrization [19], which is given as 


/+( 9 2 ) = ^E ai +) [z(q 2 ) k -(-l) k ~ K ^z(q 2 ) K ) ( 28 ) 

+ k =0 

and, 

fo(q 2 ) = - 5 - E a k )z (q 2 ) k - ( 29 ) 

^ fc =0 


Here P+.o are the Blaschke factors that take into account 
the effects of expected poles above the physical region 
but below the two body threshold f + , i.e. in the region 
(M b - M d ) 2 <q 2 < (M b + M D ) 2 , 

p+ ’° (,2)= (' - E) ■ (30) 


For / + we take the B* vector meson mass which has been 
calculated in Ref. [20] . M + = Mb * = 6.330(9)GeV. For 
the scalar form factor / 0 , there is little information theo¬ 
retically or experimentally on a 0 + bottom-charm meson. 
We take Mq to be slightly heavier than M + with large 
errors. We find that our fit results are very insensitive 
to our choice of M 0 . Even omitting the Blaschke fac¬ 
tor completely for fo(q 2 ) leads to results consistent with 


keeping it in (see test number 16 below). The poles in the 
B —)• Dlv form factors are located far above the physical 
q 2 region, for example = (M B — M D ) 2 ~ 11.6 GeV 2 
while M|„ ~ 40 GeV 2 . This implies that the form fac¬ 
tors have very small curvatures, and in fact it is very 
difficult to quantify the curvature for /o from our lattice 
data. 


Once the contributions from simple poles have been 
isolated, the power series in (28) and (29) correspond to 
smooth functions of z. One reason for prefering a power 
series in z, as opposed to one in q 2 or Ed, or even (w — 1) 
is that 1 2 : | remains very small throughout the physical 
kinematic region. For B —>• D semileptonic decays and 
our choice for to, one has 0.0 < z < 0.064. This means 
that one can go to arbitrary high powers in z k if necessary 
(in practice, with our current simulation data, going up 
to z 3 will suffice). 


The form factors of Table III are not yet in the physi¬ 
cal limit. Nevertheless, for fixed lattice spacing and pion 
mass, one can again write form factors in terms of a 
Blaschke factor mutliplying a power series in z. The ad¬ 
vantages of the z-expansion relative to an expansion in, 
for instance, powers of Ed still hold away from the phys¬ 
ical limit. What is different, however, is that expansion 
coefficients must now depend on the lattice spacing “a” 
and on “m n ” (or the light quark mass), 


(0, + ) , a(0>+) y sea \ 

a k ' y a k X ^k ( m lT m l ,°), 


(31) 


with 


c4°’ +) = x D < f’ +) (mi(phys.),m s l ea (phys.),a = 0). 

(32) 

This is the modified z-expansion first introduced in Ref. 
mm for D meson semileptonic decays, and which has 
subsequently also been employed successfully in B and 
B s meson heavy-to-light decays [SHini- The D k in (pH} 
contains all lattice artifacts and chiral logs. Specifically, 
we have, 


D k = 1 + c k x v + c k (^Sx v + Sx K ) + c k x^\og(x v ) 


+ di(am c ) 2 


dV 


2 (am c y 

2 1 Jz l 


+ e1(aE D /iry + e^aEo/ 7r) 


where 




t.K — 


Kk,v 

(4ttA ) 2 


( M \ 


asqtad ^2 
,K 


) 2 -«f Q ) 2 


(4tt /^) 2 


(33) 

(34) 

(35) 


The c k , with j = 1,2,3, and d k and e k , with i = 1,2, are 
fit parameters (we have omitted the / +i 0 label for sim¬ 
plicity) in addition to the d°’ + . In Appendix B we discuss 


what happens when the simple chiral log term in (33) is 


replaced by expressions from Hard Pion Chiral Pertur¬ 
bation Theory (HPChPT) [ 22 ] (see also test number 10 
























below). We list the priors and prior widths used in the 
chiral/continuum/kinematic extrapolation in Appendix 
C. 


We find it useful to make one more modification of the 
z-parametrization of lattice form factors. In order to ac¬ 
commodate the uncertainty coming from the truncation 
of the current matchings at 0(a s , Aqcd/M, a s /(aM)), 
we introduce new fit parameters, my and m±, with cen¬ 
tral value zero and width Smn j_, 

f\\,f± -t (l + m||)/||,(l + m_L)/j_. (36) 


The prior widths 6m\\ and Sm± correspond to our best 
estimates for higher order matching errors for Vg and 14 
respectively. With the modification of (36), our extrap¬ 
olation results coming from the modified z-expansion fit 
will then include the matching truncation errors auto¬ 
matically. To get an estimate of higher order matching 
uncertainties and fix (Smii j_> we have looked at the size 
of the known first order matching corrections. In other 
words we have gone through the correlator fits of the pre¬ 
vious section once using the fully corrected expression on 
the rhs of (12) and then a second time using just the low¬ 
est order (jjw). We find that the first order matching 
contributions have only a ~ 2% effect on fine and a ~ 4% 
effect on coarse lattices, significantly smaller than a naive 
1 x O(a) « 25 — 30% estimate. In this work we take the 
higher order uncertainties to be the same as the average 
of the full first order corrections on fine and coarse lat¬ 
tices, that is we set the prior central values and widths 
of the fit parameters mu j_ to be 0.0 ± 0.03. We have 
checked that using 0.0 ± 0.02 or 0.0 ± 0.04 everywhere, 
or 0.0 ± 0.02 for fine and 0.0 ± 0.04 for coarse lattices has 
minimal effect (see tests number 13, 14, and 15 below). 
After the modified z-expansion fits and extrapolation to 
the physical limit, these matching uncertainties for /) 
and /_l will translate into matching errors for / + and /o 
with correlations between the two form factors taken into 
account. 


In Fig. [9] we show our fit results for f + and /o plot¬ 
ted versus z. We plot both the simulation data and the 
extrapolated physical band. These are results of what 
we call our “standard extrapolation” which uses the fit 
ansatz discussed above and a z-expansion that includes 
terms through 0(z 3 ). We have carried out further tests 
of the standard extrapolation by modifying the fit ansatz 
in the following ways: 

1. stop at 0(z 2 ) in the z-expansion; 

2. stop at 0(z 4 ) in the z-expansion; 

3. add light quark mass dependence to d\ (see Eq. (30) 

of [TO]): 

4. add bottom quark mass dependence to d\ (see 
Eq. (30) of [TO]); 

5. omit (am c ) 4 term; 



z 



FIG. 9. The standard fit results with the continuum extrap¬ 
olated bands. The short horizontal bars on the upper plot 
show the fit results at non-zero lattice spacings. 


6. add ( am c ) 6 term; 

7. omit {aEjj/7 r) 4 term; 

8. add ( clEd/tt) 6 term; 

9. omit a;log(a;) term; 

10. use chiral logs from HPChPT (see Appendix B); 

11. add x 2 term; 

12. omit all Xi and xlog(x) terms; 

13. use 2% uncertainty for higher order matching con¬ 
tributions; 
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FIG. 10. Test results for /+(0) and f+(qmax) under modifica¬ 
tions of the “standard extrapolation” fit ansatz. The shaded 
horizontal bands are the standard extrapolation results. The 
x-axis labels the modifications 1-16 listed in the text. 


14. use 4% uncertainty for higher order matching con¬ 
tributions; 

15. use 2% uncertainty on fine and 4% uncertainty on 
coarse lattices for higher order matching contribu¬ 
tions; 


16. remove Blaschke factor from f 0 and f + . 

In Fig. 10 we show how results for f+(q 2 = 0) = /o(0) 
and f+(q%ta X ) are affected by these modifications. One 
sees that our extrapolations are very stable. 


V. FORM FACTOR RESULTS 


Our final results for the form factors in the physical 
limit versus q 2 are shown in Fig. 11 


f+{q 2 ) and f 0 (q 2 ) are given in Fig. 


12 


Error plots for 
We isolate the 


errors coming from different sources and also give the 
total error as a function of q 2 . The individual errors in 


Fig. 12 correspond to the following: 


• statistical 

The statistical error includes the three and two- 
point correlator fit errors and the scale errors (ri 
and ri/a). These are lattice simulation errors, and 
we have lattice data in the large q 2 re gion, from 
about 9.5 GeV 2 to 12 GeV 2 . Fig. 12 shows the 
propagation of such errors to the continuum limit 
and after extrapolation to the full q 2 range. 


• chiral extrapolation 

These are the valence and sea quark mass extrap¬ 
olation errors including effects of chiral logs. They 


FIG. 11. Continuum and chiral extrapolated /o (lower band) 
and /+ (upper band). 


conic 11 Will 




Eq. pj). 




• discretization 

Discretization errors come from the ( am c ) n and 
(■ aEjj) n terms and they constitute the dominant 
errors in our calculation. 


• kinematic 

These come from the z-expansion coefficients 
and the pole locations. As one would expect, the 
error increases as q 2 decreases. 


• matching 

Matching errors come from the mjii fit parameters 
as explained in the previous section. 


Physical meson mass input errors (0.01%) and finite size 
errors (0.1%) are not included in the plots, since they are 
too small to have any effect. 

The slope of f + (q 2 ) as one comes down from the zero 
recoil point at q 2 = qmax is a quantity that is often 
quoted when comparing different measurements of this 
form factor. In terms of the variable w = (Mg + Mg — 
q 2 )/(2MBMD) the slope parameter p 2 is given by 


Q(w) = £(1) {1 - p 2 (w - 1) + <D((w - l) 2 )} , 
where, 

G(w = w(q 2 )) = f+(q 2 ) 

1 T R 

for 

Mb 

Mb' 


(37) 


(38) 


( 39 ) 
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FIG. 12. 
/+ (upper 


Relative error components of fo (lower plot) and 
plot) for physical q 2 region. 


Another useful reference point is the value of /+(0) = 
/o(0). We find, 


/+(0) = 0.664(34). 


(43) 


In Appendix A we provide the z-expansion coefficients 
including errors and correlations for the form factors of 

Fig. [TT} 


VI. EXTRACTION OF \V cb \ 

The differential branching fraction for B —► Dlv decays 
is given by, 


dT Up\V c b\“ , m i\2\-<i 


G 2 F \Vc\ 


(44) 


(1 + ^ ' 2 #2 '“ 2 ' 


-(M 2 -M 2 ) 2 / 0 V) 


where mi is the mass of the lepton, and ijew is the 
electro-weak correction. The main goal of the present 
work is to combine experimental measurement of this 
differential branching fraction with form factors of the 
previous section to extract V c t, \. The partial branching 
fraction (the left hand side of Eq. [44 ) has been measured 
by BABAR [24] . On the right hand side, we have form 
factors from this lattice calculation, and all other factors 
are known except the target quantity |V c b|. 


In order to include the higher order electro-weak ef¬ 
fects, we apply the Sirlin factor [25], p s = 1.00662. Fur¬ 
thermore, there are final state electro-magnetic interac¬ 
tions for the neutral channel, B° —> D + lu , which we 
estimate to be a less than 0.5% effect using the signal 
yield ratio of the charged and neutral decay channels. 
Combining the two effects, we get t]ew = 1-011(5). 

We perform another modified z-expansion fit explained 
in Sec. IV together with the BABAR experiment data 
with \V cb \ as a fit parameter. We have a good fit with 
X 2 /dof = 0.88, and this is shown in Fig. 13 We get V,j } ( 
from this fit, 


\V cb \ = 0.0402(17)(13), 


(45) 


A popular way to extract p 2 is to use the Caprini- 
Lellouch-Neubert (CLN) parametrization [251 . 

2 


with, 


G(w) = ^(1) {1 - 8 p 2 z + (51 p 2 - 10 )z 2 
-(252 p 2 - 84)z 3 } , 

\/w + 1 — V2 


(40) 


(41) 


y/w + 1 + x/2 

This “z v is the same as the z-variable introduced in the 
previous section, Eq. (27), with the same to = q^ax- 

0(1) = 1.035(40). (42) 


Using Eq. (40), we extract 
p 2 = 1.119(71), 


where the first error is from the fit including all lattice 
errors and experimental statistical errors, and the second 
error is the experimental systematic error. We quote the 
experimental systematic errors as 3.3% of our fit result 
based on BABAR’s estimate of their systematic errors 
in [24]. This is equivalent to imposing 3.3% systematic 
errors on each experimental measurement bin with 100% 
correlations. 

A detailed error budget is shown in Table V. The 
dominant errors are experimental systematic, lattice dis¬ 
cretization, and operator matching errors. Thus, im¬ 
provements in both experiments and lattice calculations 
are required to obtain better precision on V c f, from our 
method. 
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FIG. 13. Form factors using both lattice and BABAR 
inputs, together with the experimental data points. 


FIG. 14. \V c b\ comparisons between inclusive and exclusive 
determinations. 


TABLE V. Error budget table for \V c t,\. The first three rows 
are from experiments, and the rest are from lattice simula¬ 
tions. 


Type 

Partial errors 

[%] 

experimental statistics 

1.55 


experimental systematic 

3.3 


meson masses 

0.01 


lattice statistics 

1.22 


chiral extrapolation 

1.14 


discretization 

2.59 


kinematic 

0.96 


matching 

2.11 


electro-weak 

0.48 


finite size effect 

0.1 


total 

5.34 



\Vcb\ has been reported from multiple lattice and non¬ 
lattice calculations. We compare the different determi¬ 
nations in Fig. [14] Our result agrees with other exclusive 
calculations, particularly with the most accurate result 
from B —> D*lv, but it is also compatible within errors 
with the inclusive determination. Since the discretization 
error is one of the dominant errors in our calculation, 
lattice errors can be reduced in the future by working on 
more ensembles with finer lattice spacings. 


leptons in the final state. BABAR has also studied decays 
involving the much heavier t lepton, B —> Dtv t , and 
measured the ratio, 


R{D) 


B{B -> Dtv t ) 
B(B -A Dlv) ’ 


(46) 


where l is either an electron or a muon. They find 


R(D)\ exp . = 0.440(58) (42), (47) 

where the first error is the statistical and the second is 
the systematic error [251 . 


Here we present a Standard Model prediction for R(D) 
based on our new form factors. Fig. 15 compares differ¬ 
ential branching fractions of Eq. (44) for B —► Dtv t and 
for B —> Dlv. Although only f+Jq 2 ) contributes to the 
Iv case, both f+{q 2 ) and fo(q 2 ) are involved in the tv t 
branching fraction. Integrating over q 2 we obtain, 


R{D)\sm = 0.300(8). (48) 

Table VI shows a detailed error budget for R(D). Fig. [16] 
gives a comparison plot for different determinations of 
R(D). All Standard Model based calculations are in good 
agreement with each other. The difference between our 
result and experiment is at the 2cr level. We note that 
we do not use any experimental results to extract R(D). 
Our result gives the most accurate pure Standard Model 
prediction to date for R(D). 


VII. THE R(D) RATIO 


VIII. SUMMARY AND FUTURE PROSPECTS 


The experimental data used in the previous section 
to extract V c j } were for semileptonic decays with light 


In this paper we have presented a new lattice QCD 
calculation of the B -a Dlv semileptonic decay form fac- 
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TABLE VI. Error budget table for R(D). 


Type 

Partial errors 

[%] 

lattice statistics 

1.24 


chiral extrapolation 

0.28 


discretization 

1.08 


kinematic 

1.61 


matching 

1.03 


finite size effect 

0.1 


total 

2.54 




FIG. 15. The differential branching fractions for B —> Dlv 
and B —> Drv 


tors f+(q 2 ) and fo(q 2 ). These were combined with ex¬ 
perimental measurements of differential branching frac¬ 
tions to extract a value for \V c b\ e xci.- Our result, given 
in Eq. 0 (and repeated in ( [45] )) is consistent with other 
recent lattice determinations using different lattice ac¬ 
tions, and provides a cross check of earlier calculations. 
We summarize these results in Fig. [14] 


The dominant error in our calculation is the discretiza¬ 
tion error, followed by higher order current matching un¬ 
certainties. The former error can be reduced by adding 
simulation data from further ensembles with finer lat¬ 
tice spacings. We are also exploring ways to improve our 
matching errors by combining simulations with NRQCD 
bottom-quarks with those employing heavier than charm 
HISQ quarks. This approach to nonperturbative match¬ 
ings of NRQCD/HISQ currents is described briefly in 
the Appendix to Ref. [TU[ . There we presented ratios 
of B s —► Klv and B s —> r] s lv form factors and explained 
how such ratios combined with a purely HISQ calculation 
in the future of B s —> rj s lv form factors would lead to 


KH 

iei 


O 

this work 

□ 

Fermilab/MILC 2015 

X 

Fermilab/MILC 2012 

V 

HQET 2010 

< 

HQET 2008 

A 

BaBar 2012 


HvH 
—<P 


R(D) 


FIG. 16. Comparisons between different determinations 
of R(D). The references for the other determinations are 
BABAR 2012 [26], HQET 2008 [27], HQET 2010 [28], Fermi- 
lab/MILC 2012 0, and Fermilab/MILC 2015 gj. 


a nonperturbative determination of the NRQCD/HISQ 
bottom-up current Z-factors. Similarly, nonperturbative 
Z-factors for bottom-charm currents used in the present 
calculation could be obtained by calculating B s —► D s lv 
forms factors once with NRQCD bottom-quarks and then 
again with heavy-HISQ bottom-quarks and then taking 
ratios. We have already completed, and are in the pro¬ 
cess of writing up, calculations of B s —> D s lv form factors 
with NRQCD bottom-quarks. Simulations with heavy- 
HISQ bottom-quarks are also underway. Hence we ex¬ 
pect to be able to significantly reduce theory errors in 
V c t, determinations from B -4 Dlv decays in the near 
future. In the meantime we hope that experimental mea¬ 
surements will also improve considerably. Only then will 
one be able to shed light on the exclusive versus inclusive 
tensions for \V c b\ via studies of B —> Dlv decays. 

In this article we also determined the ratio R(D). Our 
result is given in Eq. ([2]) (and again in (f48[). We sum¬ 
marize comparisons between Standard Model predictions 
and experiment in Fig. |16| It will be interesting to see 
whether the current ~2 a tension will develop into a true 
discrepancy between experiment and the Standard Model 
or disappear. 
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coefficient 

value 

a o°’ 

0.647(29) 


0.27(30) 

a ( 2 0) 

-0.09(2.94) 

a ( 0 +) 

0.836(33) 

a{ +) 

-2.66(52) 

4 +) 

-0.07(2.96) 


__ 

a (0) 

a o 


«2 0) °0 +) 

a l +) 

a[ +) 


a ( 0 0) 8.442e-4 -1.141e-3 -5.072e-3 

4.799e-4 3.801e-3 

5.518e-3 

a) 0) 9.255e-2 -1.087e-l 

5.390e-4 5.835e-2 

1.852e-2 

a ( 2 0) 8.652 

6.813e-3 2.504e-l 

2.402e-l 

4 +) 

1.062e-3 -7.548e-3 

-7.354e-3 

a[ +) 

2.747e-01 

-3.561e-l 

a ( 2 +) 


8.740 


simulations were carried out on facilities of the USQCD 
collaboration funded by the Office of Science of the DOE 
and at the Ohio Supercomputer Center. We thank the 
MILC collaboration for use of their gauge configurations. 


Appendix A: Reconstructing Form Factors 


We provide our z-expansion coefficients with correla¬ 
tions, so that readers can reconstruct our form factors 
for their analysis. The form factors are expressed by the 
BCL parameterization as 


um = ^ri>L + W) fc - (-1 ) k ~ K ^M 2 n 


and 


where 


fc=o 


K 




k—0 


M) = 


(Al) 

(A2) 


V t + - 9 2 - Vt+ - *0 

(A3) 

y/t+ ~q 2 + y/t+ - t 0 

{M b + M d ) 2 , 

(A4) 

drnax = { M B ~ M D f , 

(A5) 

('-*)■ 

(A6) 


For the locations of the poles, one can use M + = 6.330(9) 
GeV for / + , and M 0 = 6.420(9) GeV for / 0 to reproduce 
our form factors exactly. The coefficients, and the 

correlations are presented in Table VII. 


In the standard extrapolation of Sec. IV we used 
a generic c k x- K \og{x- K ) term to parametrize chiral loga¬ 
rithmic contributions and allowed c k to float. An al¬ 
ternate way to introduce chiral logarithms into our chi¬ 
ral/continuum extrapolations is to use expressions fixed 
by hard pion chiral perturbation theory (HPChPT) |22j . 


[!ogs]/+ = - 


2 \ r ( w ) ~ !) 


+ 1 _ 

/k ( 4 tv/tt) 

-A(xk) + A(xk) + g A(x v ) ) , (Bl) 


[ lo gs]/ 0 = - 


where 


1 T k (47r/ 7r ) 2 


(w + 1) ( r(w ) — 1) 


3 - - 1 - 

-A^) + A(x k ) + ;d(l,) 


A(x) = x log(:c), 


"( w) = 


M 2 + Ml- q 2 
2 M b M d 

log(u> + \jw 2 — 1), 


y/w 2 — 1 

k = Ad[j JAdn. 


(B2) 

(B3) 

(B4) 

(B5) 

(B6) 


We find very consistent results for the extrapolate d ph ys- 
ical form factors using either generic c k terms or (Bl) & 


(B2). This is already evident in Fig. 10 test number 10 
for /+(0) and / + (g^ ma ,). In Fig. 17 we compare the two 
approaches over the entire q 2 range. 


Appendix C: Priors and Prior Widths for the 
Chiral/Continuum/Kinematic Extrapolation 

In earlier works in mi, we split the priors for the 
modified z-expansion method into two groups: Group I 
and Group II. The Group I parameters are typical fit pa¬ 
rameters, such as quark mass dependence or z-expansion 
parameters. In this work, the Group I parameters consist 
of 


r k r k r k rt k rl k p k p k n, 

1 1 °3 ’ U 1 ! u 2 ) e l I e 2 > u k ) 


(Cl) 


where k = 0,1, and 2, and there are two sets of parame¬ 
ters for each /o and f + form factors. These parameters 
are defined in Eq. ([28]) , (l29|), and ([33]) , and the priors and 


fit results are shown in Table VIII. 


We choose priors as following. For the valence quark 
mass terms, c k , we use 0.0(1.0), since the mass terms are 
normalized by the scale, 4 . However, it is well-known 
that the sea quark mass effects are smaller than those of 
the valence quark effects, so we take 0.0(3) for the sea 
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FIG. 17. Comparison between using a generic avlog)®*-) 
term (filled blocks) and chiral logarithms from HPChPT 
(open blocks). 


quark mass terms, c§. HPChPT suggests a prior for the 
generic chiral log term, x^log(x v ), as 0.0(1). This prior 
essentially covers variations of the terms on entire kine¬ 
matic range. For more conservative error estimations, we 
take 0.0(2) as our prior for the generic chiral log term. 
We note that the prior settings with 0.0(1) and 0.0(2) 
give almost identical results. In the HISQ action, leading 
heavy quark discretization errors are 0(a s v 2 /c 2 am 2 ) and 
0(v 2 /c 2 am^). We conservatively do not take the v 2 /c 2 
terms in our power counting, so that we take 0.0(3) for 
d\ and ej priors, and 0.0(1.0) for d* and e\ priors. For 
the priors for z-expansion coefficients, we searched 

for broad enough priors that gave stable fit results, and 
we take 0.0(3.0) in this work. 


The Group II parameters are 

(a) i , <4®, aMi, (aMr tad r , (aM% sqtad ) i , 
M 0 , M+ , n , M phys , Mf ys , M phys , M p D hys , (C2) 

where i is the index for the five ensembles (i = 1 , 2,3,4, 
and 5). The Group II parameters are either from ex¬ 
periments, from other lattice simulations, or from the 
correlator fits, and used for input parameters. The prior 
settings and fit results for the Group II are shown in 
Table IX and X. 


TABLE VIII. Priors and fit results of the Group I parameters 
for the modified z-expansion fit. 


Group I 

prior [f 0 ] 

fit result [/o] 

prior [/+] 

fit result [/+] 

e? 

0.0 (1.0) 

-0.09 (18) 

0.0 (1.0) 

0.34 (20) 

cl 

0.0 (1.0) 

-0.13 (99) 

0.0 (1.0) 

-0.67 (87) 

cl 

0.0 (1.0) 

0.0 (1.0) 

0.0 (1.0) 

0.0 (1.0) 

C° 

0.00 (30) 

0.03 (28) 

0.00 (30) 

-0.10 (28) 

cl 

0.00 (30) 

0.00 (30) 

0.00 (30) 

-0.01 (30) 

cl 

0.00 (30) 

0.00 (30) 

0.00 (30) 

0.00 (30) 

cl 

0.00 (20) 

-0.10 (15) 

0.00 (20) 

0.22 (16) 

cl 

0.00 (20) 

0.006 (200) 

0.00 (20) 

0.03 (20) 

cl 

0.00 (20) 
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TABLE X. (Continued) Priors and fit results of the Group II 
parameters for the modified z-expansion fit. Parameters with 
five rows are lattice quantities for the five ensembles, Cl, C2, 
C3, FI, and F2. 


Group II 

prior 

fit result 

aM,r 

0.15990 (20) 

0.15990 (20) 


0.21110 (20) 

0.21110 (20) 


0.29310 (20) 

0.29310 (20) 


0.13460 (10) 

0.13460 (10) 


0.18730 (10) 

0.18730 (10) 

aM ^ q tad 

0.36530 (29) 

0.36530 (29) 


0.38331 (24) 

0.38331 (24) 


0.40984 (21) 

0.40984 (21) 


0.25318 (19) 

0.25318 (19) 


0.27217 (21) 

0.27217 (21) 

aM *s q tad 

0.15971 (20) 

0.15971 (20) 


0.22447 (17) 

0.22447 (17) 


0.31125 (16) 

0.31125 (16) 


0.14789 (18) 

0.14789 (18) 


0.20635 (18) 

0.20635 (18) 

Mo 

6.53 (1.00) 

6.42(43) 

M+ 

6.3300 (90) 

6.3300(90) 

ri 

0.3133 (23) 

0.3132 (23) 

M phys 

0.1373 (23) 

0.1373 (23) 

Mf ys 

0.4957 (20) 

0.4957 (20) 

A4 phy3 

5.27942 (17) 

5.27942 (17) 

A4 p ) hys 

1.86690 (40) 

1.86690 (40) 
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